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SUMMARY 

This is the second report of a set of two reports on the dynamics and control 
of slewing maneuvers of NASA Spacecraft Control Laboratory Experiment 
(SCOLE) article. In this report, the control problem of slewing maneuvers of 
SCOLE is developed in terms of an arbitrary maneuver about any given axis. The 
control system is developed for the combined problem of rigid-body slew 
maneuver and vibration suppression of flexible appendage. The control problem is 
formulated by incorporating the nonlinear equations derived in the previous report 
[1] and is expressed in terms of a two-point boundary value problem utilizing a 
quadratic type of performance index. 

The two-point boundary value problem is solved as a hierarchical control 
problem with the overall system being split in terms of two subsystems, namely 
the slewing of the entire assembly and the vibration suppression of the flexible 
antenna. The coupling variables between the two dynamical subsystems are 
identified and these two subsystems for control purposes are treated independently 
in parallel at the first level. Then the state-space trajectory of the combined prob- 
lem is optimized at the second level. 
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1. INTRODUCTION 

The primary control objective of the Spacecraft Control Laboratory Experi- 
ment (SCOLE) is to direct the RF Line-Of-Sight (LOS) of the antenna-like 
configuration towards a fixed target under the conditions of minimum time and 
limited control authority [2], This problem of directing the LOS of antenna- like 
configuration involves both the slewing maneuver of 1) e entire assembly and the 
vibration suppression of the flexible antenna-like beam. The study of ordinary 
rigid-body slew maneuvers has received considerable attention in the literature 
[3,4] due to the fact that any arbitrary large-angle slew maneuver involves 
kinematic nonlinearities. This is further complicated in the case of SCOLE by vir- 
tue of a flexible appendage deployed from the rigid space shuttle. The dynamics of 
arbitrary large-angle slew maneuvers of SCOLE model were derived in the previ- 
ous report [1] as a set of coupled equations with the rigid-body motions including 
the nonlinear kinematics and the vibratory equations of the flexible appendage. 
These nonlinear and coupled dynamical equations are m:ed in this report to study 
the slew maneuver control in terms of a hierarchical fete, back control scheme. 

The control problem of slewing maneuvers of this large flexible spacecraft is 
developed by using the two-point boundary value problem in terms of the rigid- 
body slewing and the vibration suppression of the flexible appendage as two 
separate dynamical subsystems. A decentralized optimal control scheme is utilized 
in order to solve individual boundary-value problem for each of the two subsys- 
tems by defining their state variable models and incorporating the coupling vari- 
ables between the two subsystems in these models. Also, the boundary conditions 
of the overall system are reworked in terms of bound *.ry conditions of each sub- 
system. A quadratic performance index is utilized for the overall system and is 
subsequently expressed in terms of a sum of two individual performance indices of 
the subsystems. 
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The basic algorithm for obtaining an optimal closed-loop state feedback 
scheme involves using a trajectory in terms of a vector of Lagrange multipliers as 
an initial guess at level two. This is used at level one in quasilinearization applica- 
tion. 

The two-point boundary value problem for each subsystem is solved at level 
one by using a quasilinearization technique as a trajectory optimization problem. In 
the quasilinearization procedure, starting from an ini tic. 1 guessed state trajectory, 
successive linearizations are performed of state equations in such a way that the 
final solution of the state trajectory is within an acceptable degree subject to the 
boundary conditions. The state vector definition at this level is an augmented state 
vector which includes both system states and costates. 

These optimum solutions of the subsystem trajectories are utilized at level 
two to yield the updated trajectory of the vector of Lagrange multipliers of the 
overall system to be used for quasilinearization procefs at level one. The basic 
steps of the algorithm are repeated to optimize this second level trajectory with 
respect to a prespecified error criterion to obtain an optir ial feedback law. 
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2. LIST OF SYMBOLS 

a. Vector to the point of force application on the beam 
B Damping matrix 

F 2 (.t ) Force applied at the reflector mass center 

Go (? ) Moment applied about the orbiter mass center 

I 0 Equivalent Mass moment of inertia of total assembly 

1 2 Mass moment of inertia matrix of the reflector 

J Functional used for two-point boundary value problem 

K The stiffness matrix 

L The Length of the beam 

My Effective moment applied at the reflector c.g. 

N The total number of subsystems 

Q. The generalized force vector 
qi Generalized coordinates 

JL Position vector from the orbiter mass center to the point of 

attachment 

r x x co-ordinate of the reflector mass center in the body-fixed 

frame 

r y y co-ordinate of the reflector mass center in the body-fixed 

frame 

Mi Control vector of i th system 

State vector of i th system 

Zi Vector of interconnecting variables 

X Unit vector representing the axis of rotation during the 

slew maneuver 

4> X i i th Eigenfunction corresponding to u x 

<t> yi i th Eigenfunction corresponding to iiy 

(J>yi i th Eigenf unction corresponding to uy 
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The attitude of the orbiter in the inertial frame 
Slew Angle 

The angular velocity of the orbiter in the inert al 
frame 

The angular velocity of the reflector in the inertial 
frame 

Vector of Euler parameters 
) Direc delta function 

Dual functional for two-point boundary value problem 
Vector of Lagrange multipliers 
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3. ANALYTICS 

Slew Maneuver Specification and Control Variables 

The analytics for the dynamics of SCOLE developed in reference [1] are used 
to derive the control laws for an arbitrary slew maneuver. It is assumed that the 
slew maneuver is performed by applying moments on the rigid shuttle and the 
vibration suppression is achieved by means of forces : i the flexible antenna and 
the reflector. The slew maneuver is considered to be an arbitrary maneuver about 
any given axis [1], 


The slew maneuver is defined in terms of x the axis about which the slew 
maneuver is performed. If £ is the slew angle and to is the angular velocity of the 
orbiter in the inertial frame, then the four Euler parameters can be defined as 


ci = 7isin|- 
e 2 = y 2 sin|- 
€ 3 = y 3 sin|- 
€4 = cos^- 


( 1 ) 


The four Euler parameters can be related to the angular velocity components 
of the rigid assembly as 


Cl 


Cl €4 ~ ‘C 3 £2 


0 

c 2 


C 2 C3 £4 — €1 


6>i 

C3 


€3 -€2 *1 €4 


6) 2 

C4 


C 4 -€1 -£2 -€3 


W3 


The slewing maneuver can be given in terms of the following equations [l] 


I 0 k + A 2 l = G(t)+N 2 (<*) 


(3) 
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A 2 T & + A 3 jij_ + Bg_ + K<i_ = Q(t') (4) 

where, 

G C? ) is the net moment applied about the mass center of the orbiter and is 
given by the following equation (figs. 1 & 2) 


G(f) = G 0 (r) + (r_ + a_ )xF 2 (5) 

Also, Q(t ) represents the generalized force vector which is given by the following 
equation 


L ( Qjxft ) + Qjyfr )) + Q Xl + Q yi + Q^ 

m 

.Z C Qjxp ) + Qjyft )) + q X2 + Q y2 + 


a co= 


m 

E c Q^Cr ) + Qjy t (t )) + Q xi -h Q yi + 


( 6 ) 


where, the generalized force components are given as 


L 


II 

f Fjx (z J )S(z —Zj )0 X j (z )dz 

(7) 


0 



L 


II 

£ 

6* 

f F Jy (z ,t )8(z —Zj )0yt {z )dz 

00 

w 


0 



Qjxf, t (f ) = 0 

(9) 


Here, Fj x (z ,t ) is the x component of the concentrated force applied at location j 
on the flexible antenna and F Jy is the y component of that force. 
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Also, 


<2, i O) = Z r 2*(f)'M£) 


Qy,(.t) = F 2y W)<t> yl (.L) ( 10 ) 

Qif,i it ) = M ^,(f )0^ (Z. ) 

Here, F 2 is the force applied at the reflector C. G. 

Thus, 

M^O ) = F2x r y + F 2y r x + (11) 

The location of reflector C. G. is given by coordinates (r x ,r y ) and M 2x p 
represents the external moment applied at the reflector C. G. Also, the nonlineari- 
ties N 2 can be expressed in terms of pure rigid body 1 inematic nonlinearity and 
the nonlinear coupling term between the rigid-body modes and the flexible modes. 

/Z 2 = A 4 ( 6i,0_) + A 5 ( e>,0_ ) .£. (12) 

This equation can be further simplified in terms of Euler parameters by rela- 
tionships developed in Appendix I as 

N 2 = A 6 ( + A 7 ( o>^.)i. (13) 

where j. is the Euler vector comprising all four Euler parameters. 

From equations (3) and (4) and by defining A — - A 2 Z 0 -1 A 2 + A 3 , the 
following equations are obtained 


4= A 1 


A 2 A~ 1 B<l + A 2 A~ 1 K<l + 


'fi 2 A- l A 2 T I 0 ~ l +1^ 


G(f) 


'A 2-A -1 A 2 I 0 “ 1 + Z 3 


N.i{yix£) 


(14) 
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l = -A-'Bi ~ A~ x Kg. - A~ X A 2 T I 0 ~ x G{t ) + A : A~ X A 2 T I 0 ~ X N 2 (15) 

It is assumed that control forces applied for vibrai: mi suppression has negligi- 
ble effect on rotational maneuver of the spacecraft in developing equations (14) 
and (15). Also, 1 3 represents 3x3 identity matrix in those equations. 

Subsystems and State Variable Models 

The two dynamical subsystems considered for decentralized control are the 
dynamics of the slewing of the rigidized SCOLE assembly and the vibration 
dynamics of the flexible antenna. These subsystems are represented by subscripts I 
and II respectively for subsequent analysis. 

The following are the definitions of state variables and control variables for 
subsystem I. 

*7iA«i; x J 2 A€ 2 ; x i 3 A *3; */4 A « 4 ; 

*/5 A*>i; x I6 A<o 2 ; * /7 Ao> 3 ; 


*78 AGi; Xj^jAG 2 , *710 AG 3’, 

• ♦ • 
u nA.Gi> u ^ 2 A.G 2 \ Uj$A'G$ . 

The interconnecting variables from the second subsystem to reflect the cou- 
pling between the subsystems are defined as 

Z IX A*771> ^72 A*772> z I3 A x JI3'> z IaA x IIA> ^75^-0; 

*76 AO ; z n AO ; 

^78 AO ; Z/9 A 0 ; 

^71oAO . 


The following state equations are obtained for subsystem I using these 
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definitions 


# 



X 71 





*71 

x 72 


x 72 

*73 



x 73 

*74 


0 ' D 1 0 


*74 

— 


- — I 



*75 


1 1 


x IS 

* 

= 

0 | 0 1 H 



x 16 




X I6 

J'Jl 


0 ! 0 ! 0 


X I1 

00 


X I 8 

*79 


*79 

*710 


x 710 


0 

0 


0 , 0 



I 


I 

0 I 0 


U I 1 
Ul 2 
U I 3 
~Zll 
Z I2 
Z I3 
Z I 4 


0 

H 

0 


— 2 ( x 7 1 < X I 2 > x 7 3 > x 7 4 > X I 5 > X I 6 < X I 7 ) 


Here, H = I~ 1 


a 2 a- 1 a 2 t i 0 ~ 1 + i 3 \, 


b 2 = i 0 ~ 1 a 2 a~ 1 k , 

B 3 = I 0 -1 A 2 A -1 iJ , 


( 16 ) 


D = 


x 74 — x 73 x 72 

X I 3 x 74 ~ X 71 

— x 72 X 71 x 74 

~ X I\ ~ X I2 ~ X I3 


For the second subsystem which is the flexible appendage of the entire system, 
the first two flexible modes are considered and the corresponding state variables 
and control variables are defined as 
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Minimize J (Xi»Mi»2i ) i— 1,2,...,# (18) 

w.r.t. Ui 

where Xi is the n i dimensional state vector of the ith subsystem, n* Is the 
corresponding m* dimensional control vector and z< is the r t dimensional vector of 
interconnection inputs from the other subsystem. The integer N represents the 
total number of subsystems and the scalar functional J defined by 


N 

i = 1 


V 1 

(.&(?/)) + J L t I^Cr^Cr )^(r)jdr 


(19) 


where L t | (r ) > n J - (r )^ (r ) j is the performance index at time t for i = 1,2, ..,N 

subsystems. The functional J defined in equation (19) is to be minimized subject to 
the constraints which define the subsystem dynamics, i.e. 


Xi =li 


^ ( t (f ) > z i ( t ),f 


,t 0 


( 20 ) 


Xi(O = £<0 » i = 1,2..., N . 


Also, the minimum of J must satisfy the interconnection relationship 


y Qi ( £i ) — 0 . (21) 

i= 1 

The Open-loon Hierarchical Control 

Using the method of Goal Coordination or infeasib e method [15,20], we con- 
sider another problem which is obtained by maximizing the dual function <$( X. ) 
with respect to A.(? ) ( ^ t ^ ?^ ), where 


<KD = Min 


J (x.,n 


( 22 ) 
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2L Jiil- 

subject to constraints in equations (20) and (21). Here 


*1 


Jii 


1 

--•1 

• 


• 


• 

• 

u_ = 

• 

z_ = 

• 

• 


• 


• 

£n 


i^V 




(23) 


Also, X. in equation (22) is a vector of Lagrange multipliers which is given as 


X = 


(24) 


I Aat 


~ N 

J = L 

i =1 


Pi (2j(f/)) + f Li + f X 7 ^* (*,2,,/ ) 


dt 


(25) 


j = 1,2,..., N 


Rewriting this functional J as 


~ N 

J = ZA 

i = 1 


TV 

= E 

i=l 


Pi(x i (r/)) + / jLj(x I ^^ i ,r) + X^Gf ( Xi & ,t ) 


dt 


(26) 


A 


= P i (x i (t f 



( X_i ^ ^Z_i ,t 


) + X. r G± ( x* ^ ,r ) 


dt 


(27) 


where, 
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and where Qi * ( q ,t ) has been refactored into th:. form Xf G* C x* & ,t ), i.e. 
into a form separable in the index i . 

Thus, 



(28) 


Then by the fundamental theorem of strong Lagran »e duality [22] 


min J — max <I>( X.) , i = 1,2 (29) 
Mi X. 

Thus an alternative way of optimizing J is to maximize <E> ( X. ) . 

From equation (26), for a given X.( t ) , t 0 , the functional J is 

separable into N independent minimization problems, the ith of which is given by 


Min / i =P i (x i (f/)) + J |Z* ( & Hi & ) +X. r G i (x i ^ i )J<ir (30) 

2kMi£i 

subject to 


ii = Li ). t 0 


Xi (r 0 ) 




(3D 


This leads to a two-level optimization structure where on the first level, for 
given X, the N independent minimization problems described in equations (30) and 
(31) are solved and on the second level, the X.(t)(? 0 ) trajectory is 

improved by an optimization scheme like the steepest ascent method, i.e. from 
iteration j to j + 1 
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x(ry +i = xoy +*; ±no t 0 < t/ 02) 

where 

dj = V<KAfr)) = ZGi , (33) 

<=1 

VK A) is the gradient of O ( A) > <*./ >Ois the step length and dj is the steepest 
ascent search direction. At the optimum dj -* 0 and tht appropriate Lagrange mul- 
tipier, A. > is the optimum one. 

The development of this algorithm depends on the assertion Max $( A) = 
min J and this may not be valid for all nonlinear systems. Consequently, lineari- 
zation of Q 4 , and linearized equations for may be required for constraints to be 
convex and convexity of the constraints is necessary to prove this assertion. 
Nevertheless, the method is attractive from the standpoint of simplicity and that 
the dual function is concave for this nonlinear case. This ensures that if the duality 
assertion is valid, the optimum obtained is the Global Optimum. 

On the first level, since equation (30) is to be mi 'imized subject to equation 
(31), the necessary conditions lead to a two point boundary value problem from 
which an open loop optimum control could be calculate i. However, it is desirable 
to calculate a closed loop control and for this the quasilinearization approach can 
be utilized at level one for all subsystems. Thus an iterative scheme can be set up 
whereby an initial trajectory of A^(?)* t t Q is guessed at level two and 

provided to level one. At level one the two-point boundary value problems of the 
subsystems are solved by quasilinearization. The state and control trajectories of 
all the subsystems obtained at level one are sent to level two. The test for 
optimality based on equation (33) is conducted at le- el two and if this is not 
satisfied, equation (32) is used to obtain the newX(i ) fer the next iteration. 


Subsystem Closed Loop Controllers 
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The closed loop controllers are obtained at the first level by solving the two- 
point boundary value problems of the subsystems utilizing the quasilinearization 
procedure. As noted in equations (30) and (31), the first level problem for the ith 
subsystem is 

For given X.(r ) , t 0 , 


min 




t j 

Pi |x i (f/)J + /jL i (^ji il z i ) + X 7 ’2i( J Xf^j)Jdr 


subject to 


kx — Lu (a a a ) . t 0 ^.t ^ t) 


£i ( to ) = • 


For this problem, the Hamiltonian H t can be written as 


(30) 


(31) 


Hi — Li ( x* a a ) + X r 2 i (Zali ) + Jh r Xi ( a a a ) 

For a given k, the state and costate equations become 


(34) 


a (t ) = JU ( a a a ) 


m it ) = - 


a ^ 
0a 


bLj , , 0 fJ 

6a aa “ 0a m 


(35) 

(36) 


with 


0Hi 

da 


= £L; 


0»i _ 


(37) 


It is assumed here that using the equations (36) and (37), it is possible to 
obtain the control u± and the interconnect variable vector which is an explicit 
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function of xn and Xj , i.e. 


Ui =£i ( 3k ,IU ) 


( 38 ) 


ii =sU ) 

Using these relationships for and in equations (35) and (36), the following 
equations are obtained 


i -su 

(xi,m 

). 

to <f/ 

(39) 

hi =£ 

(. 3k JU 

), 

t 0 <r r 

(40) 


with the boundary conditions 


3k ( *o ) = *io (41) 

and from the transversality conditions 


m ( tf ) = 


6 Pi 


3k ( tf ) 


a Si 


auasilinearization Procedure 


(42) 


The two-point boundary value problem of ith subsystem is given by equa- 
tions (39) and (40) subject to boundary conditions of equations (41) and (42). 
This problem is solved by quasilinearization technique as follows. 


Defining x — 


3k 


equations (39) and (40) can be rewritten as 


£(r ) = £ 


x(f ) 


(43) 
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In the quasilinearization procedure, starting from an initial guessed trajectory 
for y_=yj it ) , successive linearizations are performed of equation (43) in such a 
way that the final linear equation for y_ solves equation (43) to an acceptable 
degree subject to boundary conditions (41) and (42) which could be expressed in a 
more general form as 


X ( t 0 Y A 0 =b t> T 

(44) 

y_it f Y A f =bf r 

(45) 


where A 0 , Aj are 2 n xn matrices. 

The linearized equation of (43) about a trajectory y_ = yj it ) is obtained by 
Taylor series expansion as 

X = £(x') + ^(x')(x“X‘0 + 2P (46) 

where J (yj ) is the Jacobian of £ |y.( t ) J, t 0 ^tj , at yj and ^ represents the 
contribution of the higher order terms. Neglecting these higher order terms, the fol- 
lowing linear equation is obtained 


X = £ (x j ) + ^ • (47) 

If the initial guessed trajectory yj while satisfying equations (44), (45) and 
(47) does not satisfy equation (43), then an iterative search can be utilized to 
obtain a better linearizing trajectory by various methods discussed in references 19, 
20, and 21. This iterative search is given by noting that equation (47) can be writ- 
ten by expanding individual equations (39) and (40) by Taylor series expansion 
about a known trajectory xjit ), rj{t ), t e [ t 0 ,tf ] , and retaining terms of up to 
first order. The linearized reduced differential equations are 
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^0+D = ai (f XruJ (f)J + 


^-Cx/Cf),2I/CO) 


J§r- ( xj it ) ,ru? (r ) ) 


^0+D ( f )-^(f) 


j 

2L p+i)Cf)-r L< /(f)l (48) 


^0+i (rW (r)J + 


teL-ixSitlru* CO) 


d£j 


L ( 2^ (? ) ,n<> (r ) ) 


^0+D (f)-x/(r) 


nP +1) (f)-V(O 


(49) 


These differential equations can be rewritten as 


x^O+D = A n (t ')x i O+i\t ) + A 12 (t hi ( J +1 Kt ) + ej> (r ) (50) 

Vi (J +1) = A 2i(f )** 0 +D(r ) + a 22 (r H 0 +i )(r ) + (r ) (51) 

or, in the partitioned matrix form, 


ii^ +1) (r ) 


A n(r ) A i 2 (r ) 


XiO+iXt) 


ej. } it ) 

V'+i)(r) 


A 2i (r ) A nit ) 


i^0+i)(r ) 

+ 

Sn } Cf ) 


where the matrices 


(52) 


and 


All(r) 4-ifr * A i2 (r) 4j§-- 

A 22(OA^-, 


ej-^ A~ A n (x )xi^ {t ) - A 12 it ~)ruKt ) + & 
jLz* A~ A 2 i(r )xj i it ) — A 22 (r )t^ ^ it ) + 


(53) 


are evaluated at x± ) it ), rh J (f ) and hence are known functions of time. 

The method of complementary functions [24] can be incorporated with this 
linearization of the differential equations in the implementation of iterative search. 
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The method of complementary functions [24] can be incorporated with this 
linearization of the differential equations in the implementation of iterative search. 

An initial guess, x* 0 > 1 h° « [ t 0 , tf ], is used to evaluate matrices in equa- 

tions (53) at the beginning of the first iteration. In the next step, n sets of solu- 
tions to the 2n homogeneous differential equations 


i, 0 +1) = A u(f +1) (r ) + A tiU 0 +i)( f ) (54) 

iu 0 +1) = A 2 1 Cr )ai° +1, Cr ) + A„(t )a 0+I, Ct ) (55) 

are generated by numerical integration. For (/ +1) st iteration, these solutions are 

denoted by 2i H1 i & H2 , m H2 l m Hn • The boundary conditions 

used in generating these solutions are 


£i H1 (. ?o)=H. 21 i H Kt 0 ') = 


100 ...0 


£i H2 (. ? 0 ) = fi, ju h K r 0 )= 010. . . 0 r , 


(56) 

• • 

3U Hn it 0 ) = Q_, r 0 ) = [o 0 0 . . . 1 | r . 

Next, one particular solution at (;+l) denoted by x< p , r ^, p , is generated by 
numerically integrating equation (52) from t 0 , to tf , using the boundary condi- 
tions Xi P Qt 0 ) = Xq , rii p (f 0 ) = JQ. Then, the complete solution of equation (52) can 
be obtained by using the principle of superposition and is of the form 


i (i+1) (0 = ) + c 2 x i H Kt ) + .... + C n x i //n (r) + Xi p (.t ) (57) 

^O+D^ ) _ Cl22 .^i(; f ) + C 22 ii H Kt ) + .... + c n ru Hn (.t ) + m p (r ) (58) 

where the values of c x , c 2 c n which make 2i/ <+1 (f/ ) = U/ are to be 

determined. To find these values of Cj , c 2 c n , we let t = tf in equation 
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(58) and write it as 


21/' = 21 tu h K t/') — m Hn { */) c. + 2ii p 0/). 


C59) 


Here, c. a|ci c 2 . . . . c n 


is unknown. Solving for c. yields 


c = 


ru H Ktf ) ) 21i Hn (?/) 21/ “ 2h p • 


(60) 


It is important to note that the indicated matrix inversion in equation (60) 
has to exist in order to solve for c_. Substituting this solution of c_ into equations 
(57) and (58) gives the (y’+l) st trajectory. This completes one iteration of the 
quasilinearization algorithm and this trajectory can be further utilized to begin 
another iteration, if required. Generally, the iterative scheme is terminated by 
comparing the j th and j+1 st trajectories by calculating the norm shown in the 
following equation and comparing it with a preselected termination constatnt, p. 


Closed Loop Control 


£i 

m 


O+D 

O+D 





(61) 


In order to obtain the closed loop control, the solution of the linearized equa- 
tion (47) can be written as 


' r i 

x( t f ) = 0 ( tf ,t )x (f ) + f ,r) F (r) — J x 00 \dr (62) 

t 1 1 

where <j> is the state transition matrix of the system in equation (47). Rewriting 
equation (62) in terms of solutions of states and costates and replacing the integral 
terms by ^ (f) 



-22 - 


£i (tf ) 


011 (tf » t ) 0 12 (*/ > t ) 

Xi (?) 


£ii Cf ) 

Hi (?/ ) 


021 (?/ » * ) 022 (*/ , t ) 

Hi (?) 

+ 

£12 Cf ) 


From equations (42) and (63) 


( 63 ) 


AP 

Hi Cf/ ) = — = 021 ( ? / > (*/»*) + 022 (f/ * f )2h (f ) + Jh 2 (f ) • 


(64) 


Thus, 


Jle (f ) — 022 1 


9*< 


- 021^ (* ) ~ Pi2 (f ) 


— 022 _1 


bPi , . 

to'* 1 ™ 


022 1 I 021 Cf ) 


X t ( t ) 


(65) 


It is important to note here that 0 2 2 -1 always exists since it is a principal 
minor of the state transition matrix. 

Substituting equation (65) into equation (38) 


Hi =£i C&, t ) . 


( 66 ) 
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APPENDIX I 

The transformation that relates the orientation angles 0_ to Euler parameters _g_ 
is a nonlinear transformation. This transformation is developed for body-three 
angles representation in this appendix and similar transformations can be derived 
for other three representations, namely space-three angles, space-two angles, and 
body-two angles. 

(a) For sin0 2 1 : 

If — <0 2 <y,then 


0 2 = Sin 1 | 2 (€ 3 « 1 + € 2 c 4 ) . 

If ( cos0 jcos0 2 ) ^ 0, then 


(A.l) 


0 i = sin 1 


-2 (€ 2 € 3 - €4 ) 

cos 

sin 1 

2 ( € 3^1 + € 2 €4 ) 




If ( cos0 jCosO 2 ) < 0, then 


0 1 — 7r — sin 1 


—2 C € 2 €3 — €4 ) 

COS 

sin 1 

2( € 3 €i + € 2 € 4 ) 



If ( cos0 2 cos0 3 ) ^0, then 


0 3 = sin 1 


—2 (€i€2 ” € 3 € 4 ) 

COS 

sin 1 

2( €361 + €3^2 ) 



(A.2) 


(A.3) 


(A.4) 


If ( cos0 2 cos0 3 ) < 0, then 
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0 * = 7 t — sin 1 


—2 C e 2 — e 3 €4 ) 

COS 

sin 1 

2( e 3 €i + € 1 € 2 ) 


(A.5) 


(b) For sin0 2 = ±1, 02 is a constant. For sin0 2 =l, 0 2 = y. However, if 

sin0 2 = — 1, then 0 2 = — y . For this case, if ( sin0 !Sin0 2 sin0 3 + cos 0 3 cos0 1 ) > 0, 
then 


0 1 = sin 1 2 ( € 2 € 3 + €i e 4 ) 


(A.6) 


If ( sin0 1 sin0 2 sin0 3 + cos0 3 cos0 1 ) <0, then 


0 1 = 7 T — Sin 1 2 ( € 2 € 3 + €i €4 3 j 


(A.7) 


For this entire case, 0 3 = 0 . 
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Figure 1- Position Vectors in Inertial Frame 
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d(-L,t) 



Figure 2 - Vectors in Body-fixed Frame 



